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Abstract 

We have used the indium/copper surface ahoy to study the dynamics of surface va- 
cancies on the Cu(OOl) surface. Individual indium atoms that are embedded within 
the first layer of the crystal, are used as probes to detect the rapid diffusion of surface 
vacancies. STM measurements show that these indium atoms make multi-lattice- 
spacing jumps separated by long time intervals. Temperature dependent waiting 
time distributions show that the creation and diffusion of thermal vacancies form 
an Arrhenius type process with individual long jumps being caused by one vacancy 
only. The length of the long jumps is shown to depend on the specific location of 
the indium atom and is directly related to the lifetime of vacancies at these sites 
on the surface. This observation is used to expose the role of step edges as emitting 
and absorbing boundaries for vacancies. 

Key words: Scanning tunneling microscopy, Surface diffusion. Copper, Indium, 
Surface defects 



1 Introduction 



Over the past decade the Scanning Tunneling Microscope (STM) [1] has been 
the predominant instrument that has been employed to study atomic-scale 
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diffusion processes on surfaces. Many STM studies of surface diffusion pfie- 
nomena fiave now been performed investigating ttie mobility of steps [2-5], 
islands [6,7] and adsorbates [8]. The STM has also been the instrument of 
choice in the study of adatom diffusion and the role of adatom diffusion pro- 
cesses in crystal growth [9]. However, the study of the diffusion of naturally 
occurring adatoms and vacancies is hampered by the finite temporal resolution 
of the STM. Adatoms and vacancies both involve two energy parameters, a 
formation energy and a diffusion barrier. Typically, these are such that there is 
range of temperatures at which either species is present only in very low num- 
bers, while being extremely mobile, far too mobile to be imaged with an STM. 
Lowering the temperature to reduce the diffusion rate, so that it can be mon- 
itored with an STM, also causes the density of naturally occurring adatoms 
and vacancies to drop to extremely low values, so low that no vacancies and 
adatoms can anymore be observed. The study of adatom diffusion with Field 
Ion Microscopy (FIM) or STM [10,9] has been made possible though by the 
simple fact that adatoms can be deposited on a surface at low temperatures 
from an external evaporation source. In the case of FIM, they can also be 
produced by means of field evaporation [11,12]. No such possibility exists for 
surface vacancies and it is for this reason that the role of monatomic surface 
vacancies in the atomic-scale dynamics of single crystal metal surfaces has 
remained relatively unexposed sofar. 

In this paper we present a detailed study of the diffusion of monatomic vacan- 
cies in the first layer of Cu(OOl). For this purpose we employ indium atoms that 
are embedded in the outermost copper layer. We show that the indium atoms 
diffuse through the first layer with the assistance of surface vacancies [13,14]. 
The details of the motion of the indium contain information on the diffusion 
of the surface vacancies. The theoretical framework which we use to interpret 
our measurements is described in full in the accompanying paper, which we 
shall refer to as paper II [15]. The diffusion of vacancies leads to an unusual, 
concerted type of motion of surface atoms and causes significant mobility of 
the Cu(OOl) surface as a whole at temperatures as low as room temperature. 



2 Experimental procedures 

The Cu sample (4.8 x 4.8 x 2.0 mm^) was spark cut from a 5N-purity single 
crystal ingot. The crystal was chemically etched and then polished parallel to 
the (OOl)-plane [16]. Prior to mounting the crystal in the UHV system, we 
heated it to 1150 K for 24 hours in an Ar/H2 (20:1) atmosphere to remove 
sulphur impurities from the bulk of the crystal. After introduction into the 
vacuum chamber, the sample surface was further cleaned through several tens 
of cycles of sputtering with 600 eV Ar"*" ions and annealing to 675 K. Initially, 
after approximately every fifth cycle the surface was exposed to a few Lang- 
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Fig. 1. 662 X 404 STM image of a monatomic height step on a Cu(OOl) surface 
taken 38 minutes after the deposition of 0.03 ML of indium at room temperature. 
Embedded indium atoms show up as bright dots with a height of 0.4 A. The image 
shows a high density of embedded indium atoms near the step. (Vf = -1.158 V, It 
= 0.1 nA) 

muir of O2 at a temperature of 550 K to remove carbon contamination from 
the surface. The frequency with which the surface was exposed to oxygen was 
lowered as the preparation progressed. 

All experiments were performed with the programmable temperature STM 
constructed by Hoogeman et al. [17]. At the start of the experiments, STM 
images showed a clean, well-ordered surface with terrace widths up to 8000 
A. Small quantities of indium were deposited on the surface from a Knudsen 
cell. 



3 Vacancy-mediated surface diffusion 

Vacancies have been invoked in the past to explain the incorporation of for- 
eign atoms into a surface [18-20] or the ripening of adatom islands [21]. The 
vacancy-mediated diffusion mechanism of embedded atoms was first proposed 
for the motion of Mn atoms in Cu(OOl) during the formation of a surface al- 
loy [22,23]. Our STM investigation of the diffusion of indium atoms embedded 
within the first layer of a Cu(OOl) surface was the first to prove unambiguously 
that this motion takes place with the help of surface vacancies [13,14]. More 
recently, Pd atoms were shown to diffuse through the Cu(OOl) surface by the 
same diffusion mechanism [24] . In this section, we show how we can conclude 
that the indium atoms diffuse through the surface with the help of surface 
vacancies, without being able to observe them directly. 

The starting point of our observations is shown in figure 1. At 38 minutes after 
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the deposition of 0.03 ML of indium at room temperature, the image shows 
a high density of indium atoms in the area around a step edge. Comparison 
of the apparent height of the indium atoms observed in this image (~ 0.4 A) 
with that of indium adatoms that were observed on the Cu(l 1 17) surface (?» 
2.55 A) [25], shows that the indium atoms in figure 1 are indeed embedded 
in the first layer of the crystal. From the image it is obvious that the indium 
atoms have been incorporated in the surface through steps: the terraces are 
not populated uniformly by indium, the impurities are found only in the direct 
vicinity of the steps. This behavior is identical to what has been observed for 
indium adatoms on a Cu(l 1 17) surface [25]. Once the indium adatoms have 
reached a step and attached themselves to it, they invade the first layer on 
both sides of the step and over time slowly infiltrate the entire first layer of 
the crystal. At room temperature it takes the indium atoms typically several 
hours to spread homogeneously through the entire surface. 

The spreading of the indium atoms through the first layer implies that they 
are somehow able to diffuse, whilst remaining embedded within the surface. 
The diffusion behavior of the embedded indium atoms was studied by making 
series of images of the same area on the copper surface to form an STM movie 
of the motion [26] . 

The diffusion behavior of the indium atoms in these images is unusual in 
several respects, as is illustrated by the STM images in figure 2. 

• The jumps of the indium atoms are separated by very long time intervals. 
At room temperature these intervals can be as long as a few minutes. 

• If the indium atoms move between two images, they typically move over 
several atomic spacings. Naively one would expect them to make single, 
monatomic hops. 

• Nearby indium atoms show a strong tendency to make their jumps simul- 
taneously. If the indium atoms move independently from one another, they 
should not exhibit such concerted motion. 

To explain the unusual diffusion behavior we invoke the existence of an assist- 
ing particle. First of all, as the assisting particle is invisible in the STM-images 
we have to assume that it is too mobile to be imaged with the STM. Second, 
if the assisting particle is present in very low numbers only, this could account 
for the long waiting times. Third, a multiple encounter between the assisting 
particle and the indium atoms may lead to a jump over several atomic spac- 
ings. And finally, if the assisting particle helps to displace one indium atom, 
there is a high probability that nearby indium atoms will also be displaced, 
leading to a concerted motion of the indium atoms. It will be shown in section 
4 that also the shape of the distribution of jump lengths forms direct evidence 
for a mechanism in which the motion is only possible by virtue of an assisting 
entity. 
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Fig. 2. 140 X 70 STM images of the Cu(OOl) surfaee at room temperature illus- 
trating the diffusion of embedded indium atoms, (a) shows five embedded indium 
atoms. The right hand panel (b) shows that after 100 s the indium atoms still oc- 
cupy the same lattice sites, (c) shows the next image in which all indium atoms 
have made a multi-lattice-spacing jump. After this jump the indium atoms again 
stay at the same lattice site for another two minutes, illustrated by (d). (e-h) show 
that this pattern of long jumps separated by long time intervals repeats itself. (Vj 
= -0.580 V, It = 0.9 nA) 

We now consider the following possibilities: 

• Diffusion of embedded atoms with assistance of an adsorbed residual gas 
molecule. 

• Diffusion of embedded atoms through exchange with copper adatoms. 

• Diffusion of embedded atoms through exchange with surface vacancies. 

If the diffusion of an embedded indium atom were assisted by an adsorbed 
molecule from the residual gas in the vacuum system, the rate of long jumps 
of the indium atoms should depend on the rate at which these gas molecules 
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Fig. 3. A cross-sectional view of an exchange process of indium(bright) with a cop- 
per(dark) adatom, leading to a long jump of the indium atom, (a), A copper adatom 
arrives at the embedded indium site, either through normal hopping or exchange 
hopping(not shown), (b), The copper adatom changes places with the indium atom, 
(c) , The indium adatom now makes one or more hops over the surface before it rein- 
serts itself into the first layer. A multiple encounter between the copper adatom and 
the indium may lead to even larger dispacements. 

arc adsorbed. The diffusion mechanism could be similar to the one that was 
observed for Pt on Pt(llO) [27], where enhanced diffusion of a Ft adatom was 
enabled by adsorption of a H atom. For the diffusion of the indium atoms, such 
a mechanism can however not be active. After desorption of the residual gas 
molecule, it is no longer present on the surface and is therefore not available 
to assist other indium atoms in making long jumps. The simultaneous jumps 
of the indium atoms cannot be explained with this mechanism. Secondly, as 
we will see in section 4, the rate of long jumps depends strongly on the surface 
temperature. For the adsorption rate from the residual gas in a UHV system, 
such a dependence should not be expected. Third, the length of the long 
jumps will depend on how long the residual gas molecule resides at the indium 
atom. As the residence time of the molecule goes down exponentially with 
temperature, the jump length of the indium atom should do the same. This 
is definitely not what we observe (section 4). Diffusion of the indium atoms 
with the assistance of adsorbed gas molecules can thus be ruled out. 

This leaves the possibility of place exchange with either a copper adatom or a 
surface vacancy. For exchange with a copper adatom the diffusion mechanism 
is illustrated in figure 3. In the measurements, embedded indium atoms are 
observed to make jumps of several atomic spacings. For the adatom mecha- 
nism, this would imply that the indium adatoms reinsert themselves into the 
terrace after making at most a few hops. However, as was already seen in 
low temperature measurements [25], after deposition on the Cu(OOl) surface, 
indium adatoms do not insert themselves into the first layer. Instead they per- 
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(a) (b) 

Fig. 4. A ball-model (top view) of a diffusion event in which the passage of 
a surface vacancy leads to a multi-lattice-spacing displacement of the indium 
atom(bright). The arrow indicates the random walk pathway of the vacancy, and 
the indium- vacancy exchanges are marked with crosses to show the pathway of the 
indium between its beginning and endpoints. 

form a random walk on top of the Cu(OOl) surface and attach to steps. The 
same conclusion can be drawn from the room temperature deposition experi- 
ment, in which the indium is deposited homogeneously onto the surface. The 
adatom mechanism would lead to a homogeneous distribution of embedded 
indium atoms directly after deposition. By contrast, figure 1 clearly shows a 
high density of embedded indium atoms near a step shortly after deposition. 
The copper adatom exchange mechanism can therefore be ruled out with con- 
fidence. We identify the exchange with surface vacancies as the mechanism 
responsible for the observed diffusion of indium through the surface. Figure 4 
illustrates the vacancy-mediated diffusion mechanism. 

Having shown that the diffusion of indium in Cu(OOl) is vacancy-mediated, 
we now argue why the vacancies never show up in the STM images, by use of 
energies from embedded atom calculations [28]. Using the formation energy of 
a surface vacancy in Cu(OOl) of 517 meV, we predict a density of vacancies 
of 1- 10~^ at room temperature. Assuming an attempt frequency of 10^^ Hz, 
combined with the calculated diffusion barrier of 0.35 eV, we expect a vacancy 
jump rate of 10^ Hz at room temperature. This means that we cannot "see" 
the vacancies for two reasons. At every instant in time the probability for 
even a single vacancy to be present in the STM scan area is very low. The 
rate at which each vacancy moves is much higher than the imaging rate of 
the STM and it also exceeds the bandwidth of the STM's current preamplifier 
by several orders of magnitude. Lowering the temperature to slow the vacan- 
cies down would reduce the vacancy density even further, while increasing 
the temperature to create more vacancies would lead to an even higher va- 
cancy mobility. It is therefore only via the artificial introduction of vacancies 
above the equilibrium concentration or through the vacancy-induced diffusion 
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of tracer particles, such as the embedded indium atoms, that we can detect 
the presence and motion of the vacancies. 



4 Quantitative analysis 

The vacancy-mediated diffusion mechanism of indium atoms in a Cu(OOl) 
surface provides an unprecedented opportunity to probe the properties of 
monatomic surface vacancies. The role of the indium atom is that of the tracer 
particle. Since it can be detected with the STM it reveals precisely when a 
vacancy passed through the imaged region, and its displacement provides a 
measure for how many encounters the indium atom has had with the vacancy. 
In this section a quantitative analysis of the diffusion of embedded indium 
atoms is presented and is used to evaluate some of the fundamental energy 
parameters of surface vacancies. 

4-1 Jump length distribution 

The diffusion of indium atoms in the surface proceeds through multi-lattice- 
spacing jumps separated by long time intervals. The multi-lattice-spacing na- 
ture of the diffusion is illustrated by the jump vector distribution which is 
plotted in figure 5. The jump vector distribution shows that there is a sig- 
nificant probability for the indium atom to jump as far as five or six atomic 
spacings. In terms of the vacancy-mediated diffusion mechanism, if the va- 
cancy were making an ordinary random walk and were not influenced by the 
presence of the indium, standard random walk theory [29] gives that on aver- 
age the vacancy and the indium atom must change places as often as twenty 
to thirty times to give such a large displacement. 

In paper 11 and in ref. 14 we have derived that for the vacancy-mediated 
diffusion mechanism, one expects the shape of the length distribution of the 
long jumps to be that of a modified Bessel function of order zero. To verify 
this, the data of figure 5 were replotted in figure 6 as the radial jump length 
distribution function. Figure 6 contains six plots for the distribution at differ- 
ent temperatures. Each distribution can be fltted very well with the modifled 
Bessel function, confirming the vacancy-mediated diffusion mechanism for the 
indium atoms. Discrepancies between the expected and measured probabili- 
ties for jumps of unit length are an artifact of the automated procedure that 
was used to analyse the images [30,31]. For comparison, one of the panels 
contains the best-fit Gaussian curve that should be followed for ordinary, i.e. 
non-assisted diffusion. The only free parameter used in the fits is the recombi- 
nation probability Prec for vacancies at steps, between subsequent encounters 
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Fig. 5. The distribution of jump vectors measured from 1461 jumps that were ob- 
served in STM movies at 319.5 K. Plotted is the probabihty for jumps of an indium 
atom from its starting position to each of the nonequivalent lattice sites shown. 
Probabilities have been normalized so that the probabilities for the entire lattice 
(not just the non-equivalent sites) add up to one. In contrast to the calculations in 
paper II [15], the probability for a jump of length zero cannot be measured with 
the STM and has been put to zero. To illustrate the unusual diffusion behavior, the 
jump vector distribution for the case of simple hopping is plotted to the right. 

with the same indium atom. This probability is directly related to the average 
terrace width (see the paper II and section 4.4). The terrace widths that are 
extracted from the fitting procedure are all within a factor 2.5 of the "average" 
terrace width of 400 atomic spacings that was used in the model calculations 
of paper II. The variations in this number can be ascribed to the proximity of 
steps. The effect of steps will be discussed in more detail in section 4.4. 

4-2 Waiting time distribution 

One of the unusual aspects of the diffusion of embedded indium atoms is 
the long waiting time between consecutive jumps. The distribution of waiting 
times, expressed here in number of images, has been plotted for six different 
temperatures in figure 7. As can be seen from the figure, all measured dis- 
tributions arc purely exponential, with a time constant r that decreases with 
increasing temperature. The exponential shape of the distributions shows that 
the waiting time of an indium atom to its next long jump is governed by a 
Poisson process with rate t"^. This implies that subsequent long jumps are 
independent, which we take as proof that they are caused by different vacan- 
cies. The vacancies are created independently at random time intervals. The 
frequency with which a specific lattice site of the Cu(OOl) surface is visited by 
new vacancies is t~^. This frequency is such that our STM can easily keep up 
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Fig. 6. Radial jump length distribution for embedded indium atoms measured for six 

different temperatures. The data points have been fitted with the modified Bessel 
function of order zero that is expected for the vacancy-mediated diffusion mechanism 
[14,15]. The fits were made to the data points for jump lengths from \/2 to 6, as 
indicated with the solid part of the curve in each panel. The normalized goodness 
of fit for this range is indicated. The dashed curve in the upper left panel is the 
best-fit Gaussian curve for the same data range. 

with the diffusion process up to a temperature of 350 K. 



4-3 Temperature dependence of the jump frequency 



Simple surface diffusion processes are governed by a single diffusion barrier 
Ejj, which the difffusing species has to overcome to hop from one potential 
well to the next. To verify whether the diffusion of embedded indium atoms in 
a Cu(OOl) surface is indeed a thermally activated process, the average jump 
rate of the embedded indium atoms has been plotted as a function of l/ksT 
in figure 8. The figure shows that the diffusion of embedded indium atoms 
through the Cu(OOl) surface is indeed a thermally activated process with an 
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Fig. 7. Waiting time distributions measured for six different temperatures. All fits 

are pure exponentials The time constant r is shown in the graph for each of the 
distributions. Too high or too low count rates at short waiting times are an artifact 
of the automated analysis scheme of the images [30,31] and are ignored in the fits. 

activation energy of 699 mcV and an attempt frequency of 10^°'^ Hz. As will 
be discussed in section 6, figure 8 does not imply that the process is governed 
by a single activation energy. 



4-4 Mean square displacement 



The mean square jump length of the indium atoms is a direct measure for the 

number of encounters that an indium atom has on average with a single passing 
vacancy. From the six jump length distributions that were shown in section 
4.1, the mean square jump length has been measured. The results are shown 
in table 1. The values in table 1 exhibit a seemingly random variation that is 
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Fig. 8. Arrhenius plot of the rate of long jumps of the embedded indium atoms. The 
activation energy and attempt frequency are shown in the graph (see also fig. 11) 

much larger than the statistical error margins. As discussed in paper II [15], we 
should expect to depend only logarithmically the terrace width. Given the 
fact that all measurements were performed on terraces with a typical width 
of a few hundred Angstroms, the considerable variation in the mean square 
jump length with temperature is unexpected. 



It is at this point that the precise details of the measurements become impor- 
tant. After the deposition of the indium the room temperature measurements 
were performed and the temperature was then raised to the values shown in 
table 1. Each time after raising the temperature, the STM was allowed to 
stabilize its temperature to the point where the STM images showed no lat- 
eral thermal drift. A suitable area on a terrace was then selected to perform 
the diffusion measurements at that specific temperature. No effort was un- 
dertaken to select the same area at all six temperatures. As an example the 



T{K) 


{v?){latt.s])ac?) 


295.6 


4.37±0.11 


298.1 


4.16±0.13 


309.3 


3.61±0.04 


319.5 


7.67±0.15 


328.9 


2.82±0.02 


343.3 


5.15±0.12 



Table 1 

Measured mean square displacement as a function of temperature. 
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Fig. 9. 2495 x 1248 STM image of the Cu(OOl) surface at 343.3 K. The surface 
area that was selected to perform the diffusion measurements at tliis temperature 
is indicated by the rectangle. The trajectory to the nearest step has been drawn in 
the figure. The length of this trajectory is shown in table 2 (Vt = -1.158 V, It = 0.1 
nA). 
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Table 2 

Measured nearest-step-distances as a function of temperature. All distances have 

been measured in Cu(OOl) atomic spacings. The two values marked with a * are 
estimated values: the overview scan at 328.9 K showed only a part of the relevant 
surroundings of the area that was used to obtain the jump data. 

region that was selected for the 343.3 K measurements is shown in figure 9. 
From each of the overview scans that were made prior to starting the diffusion 
measurements, the distance to the steps surrounding the measurement site 
was quantified by measuring both the logarithm of the minimum distance to a 
step as well as the polar average of the logarithm of the distance to the steps. 
The measured values are shown in table 2. 

Both ways to quantify the proximity of steps to the measurement site show 
that there is a very clear trend for longer jumps when the steps are further 
away. This trend is even more clear from figure 10, where the mean square 
jump length is plotted versus the average logarithm of the distance of the 
measurement site to a step. In this graph the 328.9 K measurement has not 
been used as the measurement area was partly outside the area that was 
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<log[nearest step distance]> 

Fig. 10. Plot of the mean square jump length versus the average logarithm of the 
distance of the measurement area to the steps. The graph shows that measurements 
that are made further away from a step yield larger jump lengths. All distances were 
measured in atomic spacings. 

imaged in the overview scan and the distance to the surrounding steps could 
not be properly determined. 

Because all measurements were performed in thermodynamic equilibrium, the 
density of vacancies throughout the terraces is constant. The observation that 
the jump length of the indium atoms depends on the specific geometry in 
which the measurements were performed therefore implies that the rate of 
long jumps is also position dependent. This in turn means that since all mea- 
surements were performed at different areas, the temperature dependence of 
the diffusion of embedded indium atoms can only be properly measured by 
plotting the diffusion coefficient D of the indium atoms versus l/ksT as op- 
posed to the jump rate versus l/ksT which was plotted in figure 8. This has 
been done in figure 11. We observe that the scatter that was present in figure 
8 is significantly reduced and a much more accurate value of 717 ± 30 meV 
for the activation energy is obtained. 



5 A vacancy attachment bcirrier ? 

Indium atoms that are incorporated into the first layer of the Cu(OOl) surface 
at room temperature via a step appear to have a slight preference for the upper 
terrace with respect to the lower terrace. This effect can be seen from figure 
1 where approximately 60 % of the indium atoms have entered the upper 
terrace and only 40 % have gone into the lower terrace. The effect is even 
more dramatic if the incorporation of the indium takes place at temperatures 
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Fig. 11. Arrhenius plot of the diffusion coefficient D of embedded indium atoms. 
The coefficient D was calculated at each temperature by multiplying the rate of 
long jumps (figure 8) with the mean-square jump length (table 1). 




Fig. 12. 830 X 415 STM image of indium being incorporated into the upper 
terrace on the right hand side of a step on the Cu(OOl) surface. The image was 
taken at a temperature of 280.4 K during a slow ramp of the temperature from 130 
K, the temperature where 0.015 ML of indium was deposited on clean Cu(OOl). The 
indium first started to be incorporated into the upper terrace at 267 K and shows 
up as a light band on the right hand side of the step. (Vj = -0.672 V, Ij = 0.1 nA) 

below room temperature. Figure 12 shows an image that was obtained at 
280.4 K, after indium had been deposited on the surface at 130 K and was 
first observed to be incorporated into the surface at a temperature of 270 
K. The image shows that practically all the indium atoms enter the upper 
terrace at this temperature. This observation may imply the existence of a 
step edge barrier, similar to the Schwoebel barrier for adatoms [32-34], but 
in this case for the attachment of surface vacancies at the lower side of a 
step. Using the ratio of the number of indium atoms that is incorporated 
into the upper and lower terraces at 270 K and 300 K, we derive a first 
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Fig. 13. One-dimensional energy landscape for a vacancy in the Cu(OOl) surface. 
For this one-dimensional situation the vacancy is formed at the step on the left and 
approaches the indium from there. It is only able to displace the indium atom to the 
left by one atomic spacing. Ep is the formation energy of a Cu(001)-vacancy, Ed 
is the diffusion barrier for a Cu(001)-vacancy, Ej^ is the vacancy-indium exchange 
barrier and Ejn is the binding energy for an indium- vacancy pair. 

estimate of this barrier of 0.8 eV. The difference in incorporation between 
the upper and lower terrace has already been discussed in the context of the 
Mn/Cu(001) surface alloy [35,36,23] and can in principle be attributed to 
the difference in the incorporation processes for the upper and lower terrace. 
However, the high value of 0.8 eV suggests that other factors play a role in 
the preferred incorporation into the upper terrace. One such factor could be 
the large amount of indium that is present at the step during incorporation. 



6 Interpretation of the activation energy 

By definition, the rate at which the indium atom is displaced by a surface 
vacancy is the product of the vacancy density at the site next to the indium 
atom times the rate at which vacancies exchange with the indium atom. The 
activation energy obtained from the temperature dependence of the total dis- 
placement rate will yield the sum of the vacancy formation energy and the 
vacancy diffusion barrier. When the measurements are performed with a finite 
temporal resolution and if there is an interaction present between the vacancy 
and the indium atom, this simple picture changes. 

The energy landscape for the vacancy is sketched in figure 13. We now consider 
the situation where the diffusion is observed with a finite time resolution. We 
assume that the time resolution of the STM-measurements is sufficient to 
resolve the effect of the entire random walk of one individual vacancy from 
that of the next, which was already shown to be the case in section 4.2. The 
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Fig. 14. The random walk starting situation, with the vacancy directly next to the 
indium atom. The probabilities x, y and p^ec are introduced in the text. 

starting situation is shown in figure 14: an indium atom is embedded in the 
origin of a square lattice with a vacancy sitting next to it at (1,0). The only 
diffusion barrier that has been modified by the indium atom is the one which 
is associated with the vacancy-indium exchange. This is equivalent to putting 
Ein = in figure 13. 

We now define the following additional quantities: 

X the probability that in the first step the vacancy imme- 
diately changes places with the indium. 

y the probability that the vacancy docs not immediately 
change places with the indium but that it returns to any 
of the four sites neighbouring the indium. 

Prec the recombination probability of the vacancy, i.e. the 
probability that the vacancy recombines at a step with- 
out ever changing places with the indium. This proba- 
bility is determined by the distances to the nearby steps. 

Note that x -\- y + Prec 7^ 1, because, as is indicated in figure 14, there is the 
possibility that the vacancy returns one or more times to the sites neighbour- 
ing the indium atom and then recombines, without ever exchanging with the 
indium. Only if one were to define the recombination probability as p^g^, the 
probability that the vacancy recombines without first returning to any of the 
sites neighbouring the indium atom, would x + y+p'^.^^ — 1. The rate at which 
jumps of the indium are registered by the STM is determined by the rate at 
which total displacements with a length of at least one atomic spacing occur. 
Given the starting situation of figure 14, the probability that the indium is 
displaced by the vacancy is given by 

17 °° X 

Pdis = 1 - Prec ^ X + y ■ X + y^ ■ X + ... ^ Y^y^ .x= —— (l) 

N=Q y 



Here the possibility that multiple displacements of the indium may add up to 
a zero net displacement is ignored[^ The precise value of a:, the probability to 
make a successful exchange, is determined by the energy difference between 
the modified and unmodified vacancy exchange barrier, Eu — E^. 



X = x{T) 



Z/qC 



z/qc ''b'^ + 3z/oe '=sT 1 + 3e 



Upon examining the random walk of the vacancy when it has not exchanged 
with the indium but has returned to any of the sites neighbouring the indium, 
we see that only the first step of this random walk pathway has a temperature 
dependent probability associated with it and y can therefore be split up in a 
temperature dependent and a temperature independent part 

y = y{T) = {I - x{T))C (3) 



where < C < 1 is a constant which is determined by the geometry of the 
lattice, and especially by the distribution of distances between the starting 
position of the vacancy and the absorbing boundaries (steps). It can be eval- 
uated numerically. Substituting equation (3) in (1), the probability to have at 
least one displacement is equal to 



Pdis 



1 - yiT) 

x{T) 



l-{l-x{T))-C 
x{T) 



(1 - C) + Cx{T) 

^ 7—r- (4) 



D. 

T 



1 + 3(1-^)6 

From equation (4) it is clear that the final rate of long jumps will contain 
exponential terms not only in the numerator, but also in the denominator. 
The rate of long jumps should not show normal thermally activated Arrhenius 
behaviour. 



The observed rate of long jumps is equal to the equilibrium rate at which 
vacancies exchange with the indium atom, divided by the average number 

^ This is justified by numerical calculations presented in paper II where we address 
this problem for the case of indium. Our calculations show that the fraction of jumps 
that have a zero net length changes only marginally with temperature. 
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of elementary displacements caused by a single vacancy, given that the va- 
cancy has displaced the indium atom at least once. This average number of 
displacements (n) is given by 



(n) = 1 + • Prec + 1 • (1 - Prec)Prec + 2 • (1 - Precfprec + ■■■ 
1 



Prec 

1 / 1 



-C] (5) 



1-C \l-x{T) 
Using equation (5), the observed rate of long jumps is equal to 



^LJ = 1^0 (6) 



1 + 



e "B^ 
C 



where 
C = 3(1 - C) 

In the case that indium is "identical" to copper, so that A = 0, ( n ) becomes 
3{i-c) ^^^^ long jumps reduces to 

e ''bT y 
i -t- ^, C' 



where v is the jump rate of atoms in the clean Cu(OOl) surface. The denomi- 
nator corresponds to the number of exchanges contributing to a long jump. 

In the extremely repulsive case 3> — 1 and the average number of dis- 
placements (n) becomes approximately equal to 1. Equation (6) reduces to 

^LJ = i^oe ''bT = VqB ''bt (8) 



In the extremely attractive case <^ 1 and the average number of jumps 

A 

(n) is ■ \e ''bt . Equation (6) now reduces to 

VLJ = 3i/o(l - C)e ^bt (9) 



19 



We immediately see the fundamental difference between attraction and re- 
pulsion. The apparent activation energy for strong attraction is identical to 
that for the copper surface itself, while it is larger in the case of repulsion. 
This asymmetry between attraction and repulsion is caused by the fact that 
for moderate to strong attraction, the probability x rapidly approaches unity, 
meaning that the arrival of a vacancy next to the indium is almost guaranteed 
to cause a long jump. In contrast, in the case of repulsion, the probability 
X scales with the Boltzmann-factor containing the exchange barrier of the 
vacancy and the embedded atom. 

In general, the Arrhenius plot of the log of the rate of long jumps versus 
is nonlinear: 

In {vLj) = In (i/o) - (^F + + A) -i- - In I 1 + 1 (10) 



The last term on the right describes the departure from ideal Arrhenius be- 
haviour. With this term expanded to second order in -j^, the rate of long 
jumps can be rewritten as 



In (i/Lj) = ln I UoT^TT, \ + [-Ef-Ed+[ rr^, - 1 ) A 




1 + C" ) ) ksT 



2{l + C'Y{kBT) 



For practical values of we are in the limit <^ — 1 and the measured 
activation energy is equal to the sum of the vacancy formation energy, Ep^ 
and the vacancy diffusion barrier, Ejj. 

We mention again that the fact that multiple displacements of the indium 
may actually add up to zero is ignored. The numerical calculations in paper 
II [15] show that for indium approximately 8.5 % of all long jumps have a net 
displacement of zero atomic spacings. This fraction varies over the investigated 
temperature range by an insignificant amount and will therefore only result 
in a small modification of the prefactor in the temperature dependence. 



7 Discussion 



The jump length distributions mentioned in this work all have the shape that 
is expected for the particle-assisted diffusion mechanism. All distributions can 
be accurately fitted with the modified Bessel function, reconfirming that the 
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diffusion of indium through the first layer takes place with the aid of some 
particle, which we deduced to be a vacancy by looking in detail at the incorpo- 
ration process of the indium atoms. Having used the jump length distributions 
to validate the model described in the accompanying paper, we can now use 
this model to calculate the diffusion coefficient of copper atoms in a clean 
Cu(OOl) surface. This involves no more than taking the interaction between 
the vacancy and the indium atom out of the model such that the vacancy 
performs an unbiased random walk. At room temperature the average jump 
length of copper atoms in the clean Cu(OOl) surface is calculated to be 1.6 
atomic spacings. Using the rate of long jumps of the indium atoms that was 
measured at room temperature we find that the diffusion coefficient of surface 
atoms in the Cu(OOl) surface as a result of vacancy mediated diffusion is 0.42 
A^-s~^. Using the attempt frequency that was obtained from figure 8, we find 
that at room temperature every site on the surface is visited by a new vacancy 
on average once every 32 s. 

The sum of the Cu(OOl) surface vacancy formation and migration energy is 
equal to 717 ± 30 meV. The sum of the EAM calculated diffusion barrier (0.35 
eV) and the formation energy (0.52 eV) amounts to 0.86 eV. The fact that 
this result is too large is not entirely unexpected as all diffusion barriers that 
were calculated with EAM were too large when compared with experiments. 
However, the value of 0.35 eV was already obtained by dividing the true EAM 
value by a factor 1.7 [28]. Using the measured sum and the calculated EAM 
formation energy of 0.52 eV, we find a lower value of 0.20 eV for the diffusion 
barrier of a Cu(OOl) surface vacancy. Using the vacancy formation energy 
of 0.485 eV from [24] that was obtained through first-principles calculations 
we obtain a vacancy diffusion barrier of 0.232 eV. With such a low diffusion 
barrier, we expect a room temperature vacancy diffusion rate of 10^ Hz. 

The interpretation of the activation energy as the sum of the vacancy diffusion 
barrier and the vacancy formation energy is supported by measurements of 
vacancy-mediated diffusion of Pd in Cu(OOl) [24], which reveal a value of 
0.88 ± 0.03 eV for the sum of the vacancy formation energy and the Pd- 
vacancy exchange barrier. The higher activation energy is a direct consequence 
of the fact that vacancies and embedded palladium atoms repel (see section 
6). To be able to experimentally separate out the two energy parameters, 
the vacancy formation energy and the vacancy diffusion barrier, independent 
measurements are needed of one of the two energies. One possibility is to follow 
artificially created surface vacancies at low temperature [39,37,38]. Such a 
measurement would yield the diffusion barrier for a monatomic surface vacancy 
which could then be combined with the 717 meV that we measured to obtain 
the vacancy formation energy. 

The variation in mean square jump length that was observed with the distance 
to the nearest step confirms the idea that surface vacancies are formed and an- 
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nihilated at steps. The creation of a surface vacancy through the expulsion of 
an atom from a section of flat Cu(OOl) terrace has been predicted to cost 984 
meV [40] . Because of this high energy for the creation of an adatom- vacancy 
pair, at most temperatures surface vacancies are formed almost exclusively 
at steps. This means that the steps section the surface into separated areas 
through which the vacancies are allowed to diffuse. The vacancies are anni- 
hilated when they reach a step. A vacancy which exchanges with an indium 
atom near a step has a relatively high probability to recombine at the step 
as opposed to making another exchange with the indium atom. Indium atoms 
near a step will therefore on average make shorter jumps than indium atoms 
that are far away from a step. Given a constant density of surface vacancies 
throughout the terrace (the measurements were performed in thermal equilib- 
rium), this implies that the rate of long jumps near steps will be higher than in 
the middle of a terrace. Hence, the indium will perform relatively many jumps 
of a short length near a step, whereas in the middle of a wide terrace, indium 
atoms will tend to jump less frequently, but given the increased lifetime of the 
vacancy, they will jump over longer distances. This observation opens the way 
to detailed studies of vacancy creation and annihilation near steps and kinks. 

The existence of the vacancy attachment barrier that was proposed in section 
5 needs to be investigated in more detail. The influence of the considerable 
amount of indium that was present at the step makes it hard to draw detailed 
conclusions about this barrier. Further experiments investigating the existence 
and magnitude of this barrier at a clean step should be performed with low 
densities of embedded guest atoms, e.g. in adatom or vacancy islands. The 
lifetime of a surface vacancy in such islands will be affected strongly by the 
presence of this barrier. If present, the barrier will lead to significant differences 
in the mean square jump length and frequency of long jumps of indium atoms 
embedded in equal-size adatom and vacancy islands. 



8 Conclusions 

Surface vacancies were shown to be responsible for the motion of embedded 
indium atoms in the Cu(OOl) surface. The density of surface vacancies at room 
temperature is extremely low, but vacancies are able to diffuse through the 
surface at an extremely high rate. In the STM-measurements the rapid diffu- 
sion of these vacancies leads to simultaneous long jumps of embedded indium 
atoms. Measurements of the jump length and jump rate were performed at six 
different temperatures. These measurements show that the vacancy-mediated 
diffusion process can be accurately described with the model that is presented 
in paper II [15] and ref. 14, provided that the interaction between the em- 
bedded indium and the surface vacancies is properly taken into account. The 
role of steps as sources and sinks for surface vacancies has been confirmed by 
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measuring the position dependent jump rate and jump length of the indium 
atoms. Near a step the jump rate of the indium atoms is increased, but at 
the same time, given the shorter hfetime of vacancies near a step, the jump 
length is decreased. The sum of the vacancy formation and migration energy 
was measured to be 717 ± 30 meV. 

The most far-reaching conclusion concerns the mobility of the Cu(OOl) surface 
itself. The diffusing vacancy leaves behind a trail of displaced surface copper 
atoms and thereby induces a significant mobility within the first layer of the 
surface, already at room temperature [13]. The mobility of atoms in the clean 
Cu(OOl) surface was evaluated by applying the model of ref. 14 and paper 
11 [15] with the interaction between the tracer particle and the vacancy ignored. 
Through this procedure the diffusion coefficient of Cu(OOl) surface atoms was 
measured to be 0.42 A^-s~^ at room temperature. 
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